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ABSTRACT 



Radio astronomical measurements of extended emission require knowledge of the 
^ ■ beam shape and response because the measurements need correction for quantities 

2^ ■ such as beam efficiency and beamwidth. We describe a scheme that characterizes the 

^T) . main beam and sidelobe in all Stokes parameters employing parameters that allow 

t^ I reconstruction of the complete beam patterns and, also, afford an easy way to see how 

the beam changes with azimuth, zenith angle, and time. For the main beam in Stokes 
/ the parameters include the beam width, ellipticity and its orientation, coma and 
1-^ ■ its orientation, the point-source gain, the integrated gain (or, equivalently, the main 

I ' beam efficiency); for the other Stokes parameters the beam parameters include beam 

^ ■ squint and beam squash. For the first sidelobe ring in Stokes / the parameters include 

c/3 . an 8-term Fourier series describing the height, radius, and radial width; for the other 

Stokes parameters they include only the sidelobe's fractional polarization. 

We illustrate the technique by applying it to the Arecibo telescope. The main 
^ ' beam width is smaller and the sidelobe levels higher than for a uniformly-illuminated 

C^ ' aperture of the same effective area. These effects are modeled modestly well by a 

blocked aperture, with the blocked area equal to about 10% of the effective area (this 
corresponds to 5% physical blockage). In polarized emission, the effects of beam 
squint (difference in pointing direction between orthogonal polarizations) and squash 
(difference in beamwidth between orthogonal polarizations) do not correspond to 
theoretical expectation and are higher than expected; these effects are almost certainly 
caused by the blockage. The first sidelobe is highly polarized because of blockage. 
These polarization effects lead to severe contamination of maps of polarized emission 
by spatial derivatives in brightness temperature. 



Subject headings: polarization — instrumentation: polarimeters — techniques: 
polarimetric — techniques: miscellaneous 



1. INTRODUCTION 

The AreciboF] radio telescope is unique because the feed moves with respect to a reflector 
that is fixed with respect to the ground. Thus, the properties of the telescope beam can change in 
a complicated way as a source is tracked because some components of the beam rotate with the 
parallactic angle and some don't. For example, the beam is significantly elliptical, and the ellipse 
rotates with parallactic angle; in contrast, the effects of reflector surface imperfections are fixed on 
the ground. Disentangling these complicated behaviors is best done by parameterizing the beam: 
some parameters change in recognizable and physically reasonable ways. 

Our efforts were directed towards characterizing all of Arecibo's receivers with the partial 
intent of establishing a historical record so that the beam patterns could be monitored as future 
improvements in the telescope structure occur. However, our efforts are also useful for the process 
of interpreting observational results, and can easily be applied to other telescopes. This paper 
is a condensed version of Arecibo Technical and Operations Memo 2000-03, which covers all of 
Arecibo's receivers and with more detail. 

We devised a simple and fast observing pattern that provides nearly complete information on 
the main beam and first sidelobe in all four Stokes parameters. We use the digital crosscorrelation 
technique to derive the Stokes parameters (Heiles 2001). The pattern consists of four continuously 
sampled scans oriented every 45° with respect to {dz, za) yielding an 8-pointed star-shaped pattern 
that extends out to three nominal HPBW (Half Power Beamwidths) on each side of the nominal 
pointing position. 

These data allow parameterization of the main beam and first sidelobe as described in the 
following sections. We begin in §|2| by describing the parameterization scheme for Stokes /. §^ 
describes representative Stokes / results for the main beam parameters. ^ compares those results 
with a theoretical model of a uniformly-illuminated blocked circular aperture. §^ presents Stokes 
/ results on the first sidelobe. §^ defines and describes the parameters beam squint and squash 
for the polarized Stokes parameters, presents the polarization of the first sidelobe, and calculates 
the levels of contamination caused by spatial variations in brightness temperature on maps of 
polarized emission. §0 presents greyscale images of the beam in all four Stokes parameters. 



PARAMETERIZATION OF THE STOKES I BEAM 



^The Arecibo Observatory is part of the National Astronomy and Ionosphere Center, which is operated by Cornell 
University under a cooperative agreement with the National Science Foundation 



2.1. Five Parameters that Describe the Main Beam Shape 

We follow Rohlfs and Wilson's (2000) definition of the normalized main-beam power pattern, 
P{9,(j))n = p ■ We represent it with five parameters of a two-dimensional Gaussian defining a 
power pattern with ellipticity and coma. The parameters are the average beamwidth Bq, beam 
ellipticity 0i with an orientation ^beami a-nd coma Ocoma with orientation (/'coma- 

To parameterize the elliptical beam, we take the ratio of major to minor axes to be q°J^q^ 
with the major axis oriented at position angle (/'beam! and the beamwidth Q to vary as 

G = Oo + ©1 cos 2((/) - ^boam) • (1) 

where (p is the position angle in the {az,za) coordinate system, defined to be zero along the 
positive azimuth axis and increase towards positive zenith angle (see Figure |l|). 

We parameterize the coma lobe with a function chosen for analytical convenience; it is not 
based on a proper derivation. The parameters are the coma strength Ocoma and its position angle 
'/'coma- We assume that the effect of coma depends only on the projection of 9 along the coma 
direction; this distance is 



6'coma = COs{(j) - (/'coma) - (2) 

The normalized power distribution including parameterized ellipticity and coma is then: 



Pn = exp 



^2 l-min{Qcoma%r,0.75} 



02 



(3) 



where 9 is the great-circle angular distance of the position from the true beam center, equal to 
Wa^^g-gg^. + ^^flofj-sef '^^^ term minjocoma q"" ,0.75} is included to prevent the coma term from 
unduly distorting the beam far from the center. Note that the beamwidth Q related to the usual 
half-power beamwidth by HPBW = 2(ln 2)^20 ^ 1.6659. 

Figure || illustrates these quantities with a main beam having severe coma and exaggerated 
ellipticity. In this example, Ocoma = 0.2, (/'coma = 22.5°, Qq = 3.4 arcmin, 0i = 1 arcmin, and 
(/'hpam = -67.5°. 



2.2. Relating the Beam to the Observation Pattern 

When observing the beam shape we must point the telescope to an assumed beam center. 
The four scans of the observing pattern intersect this assumed beam center, as shown by the 



dash-dot lines on Figure |l[ The difference between the assumed and actual centers is the pointing 
offset. Thus, in equation ^, 

aZoSsct = azohs + Saz (4a) 

^aoffset = zaohs + Sza (4b) 

where (a-Zobsi -^aobs) are the "observed" offsets, that is the ones measured with respect to the 
assumed center. 



2.3. The First Sidelobe in Stokes I 

Detailed maps of Arecibo's first sidelobe show some azimuthal structure (i.e., with (p on 
Figure ||) (Heiles 2000). Our observing pattern provides four radial cuts through the beam. We 
least-squares fit each radial cut with three Gaussians, two weak ones at the beginning and end 
for the sidelobes and a strong one for the main beam. This yields values of the sidelobe height, 
center, and width at 8 (/)'s. For each one, e.g. the height, we Fourier transform the 8 values and 
derive 8 complex Fourier coefficients, only 5 of which are independent because the inputs are real 
numbers. We reconstruct the sidelobe using A^ Fourier coefficients, with N < 8. Using N < 8 
is not unreasonable because the measurement and Gaussian fitting processes are hardly perfect; 
it is equivalent to performing a least-squares fit on the 8 (/)'s and deriving fewer coefficients. In 
practice, we reconstruct the sidelobe with A^ = 6, which is just sufficient to reveal the blockage 
caused by the equilateral triangle feed support structure, which is a major structural component. 

This procedure is complicated in practice. In some cuts the sidelobe is absent or is not 
well-fit. We test each derived sidelobe Gaussian fit and accept it only if its width is smaller than 
the nominal HPBW&nd greater than 0.3 times the nominal HPBW . We also require the fractional 
uncertainty in the derived intensity to be less than 0.45. In cases having no acceptable fit, we 
take the height equal to zero, and take the center and width equal to the average of all the other 
acceptable centers and widths for that particular group of cuts. 



2.4. Effective Area, On- Axis Gain, and Beam Efficiencies 

Again we follow Rohlfs and Wilson (2000) and define the directive gain G as the gain relative 
to an isotropic radiator for a single polarization. Here we are concerned with treating polarization, 
for which purpose we must explicitly define the unsubscripted directive gain to apply to ^'^'^ i 
below we will subscript it with a polarization label. Conservation of energy requires that G satisfy 
the relationship 

Gdn = Air . (5) 

whole sky 



The on-axis point-source directive gain is directly related to the effective area ^eflf by 

and in practical units of source flux in Jy, area in m^ we have 

KprJy = W'^^^ (7) 



or 



^eff = 2 X W'''kBKperJy (8) 



so, eliminating A^s 



_ 47r(2 X lO^^kBKpsrJy) 

3. 47 X 10^ KperJy 

A2 ^ ' 

where Acm is the wavelength in cm. Arecibo's theoretical gain, without blockage, is KperJy ~ 11, 
which is equivalent to a uniformly illuminated circular aperture with radius 94 m. The normalized 
power pattern P„, which we used above in equation |3|, is just 

Pn = yr— = 2-89 X lO-'^r^ , (11) 

<-^max KperJy 

so 

Pndn = 3.62 X IQ-^-^SV • (12) 

whole sky KprJy 



J J V 



The "efficiency" of a portion of the telescope's beam is defined as the integral of its directive 
gain divided its integral over the whole sky. Thus, using units of arcmin^ for Q and cm for A, 

Jyupt Jij r p 

^beam = /.o4 -^ / / i->jasiarcmin2 ■ (l"j) 

'^cm -^ •' beam 

Normally we consider the main beam efhciency, but we can also use the same terminology for the 
sidelobes. The main beam and first sidelobe integrals are, in principle, expressible analytically by 
our parameterization scheme; in practice, however, we compute them numerically. 



3. REPRESENTATIVE OBSERVED PARAMETER VALUES 

Table |l] provides representative observed values for the above-defined parameters for two 
receivers at Arecibo Observatory, measured during September 2000. Both are point-source feeds 
used with the new Gregorian three-mirror system. The point-source gain at 430 MHz is close 



to that expected for a perfect telescope, while in for the L-band Wide receiver (LBW) the gain 
decreases with frequency. This is almost certainly a result of irregularities in the main reflector; a 
program of adjustments is planned for early 2001. 



3.1. Main Beam Ellipticity and Coma 

Arecibo's beam is elliptical by design; (/)beam should be 90°, meaning that the za beam is 
larger than the az beam. In practice we find that (j)beam exhibits systematic departures from 90°, 
ranging from ~ 80° — > 110° and does not get noticeably worse with frequency. We believe this is a 
result of blockage. 

In September 2000 Arecibo had a noticeable coma lobe. It gets worse with frequency, with 
strength Ocoma rising from about 0.015 to 0.15 from 430 to 4500 MHz (at za ~ 10°). The direction 
'/'coma exhibits no clear systematic departure from being constant as a function of [az, za) or 
frequency, scattering about (j) ~ 25°. This suggests that it is tied to the azimuth arm and therefore 
probably a result of subreflector misalignment. 



4. COMPARISON WITH A THEORETICAL MODEL: A 
UNIFORMLY-ILLUMINATED CIRCULAR APERTURE HAVING 

BLOCKAGE 

It is instructive to compare Arecibo's beam pattern with a theoretical model. The most 
convenient is a uniformly-illuminated circular aperture. We will then consider a circular aperture 
with a central circular blockage. Even though Arecibo's beam is elliptical, it's close enough to 
circular that the model makes physical sense. 



4.1. A Uniformly-Illuminated Unblocked Circular Aperture 

In this section we will define and evaluate the theoretical parameters for an unblocked 
aperture, with which we will compare our corresponding observed parameters. This comparison 
leads to a profusion of subscripts, and to avoid confusion we will first define those subscripts. The 
subscript "?7" denotes the theoretical result for a uniform illumination pattern. There are two 



forms of such results. One is obtained directly from the theory beginning with equation |lj below, 
and this is denoted with the additional subscript "/". The other is an approximate theoretical 
result obtained by fitting with Gaussians (as we do the observational data) the aforementioned 
direct theoretical result (subscripted with "I"); this approximation is the one we compare with 
the data, and has no additional subscript. In addition, the subscripts ^''MB" and "F5" refer to 
the main beam and first sidelobe, respectively. 



where 



The electric field pattern of an unblocked aperture is 

F(sm9) = 2:^^, (14) 

u 

n=^ and 0, = A/Z? 

(Rohlfs and Wilson 2000). The field for a blocked aperture is a suitably- weighted arithmetic 
combination of different diameters. In principle, the best way to fit this model to the data is 
to compare the observed and theoretical beam shapes. However, this is not so easy because of 
Arecibo's elliptical beam. Rather, we examine five parameters. Of these, two are robust in the 
sense that they don't depend on assuming a perfectly smooth refiecting surface for the bowl and, 
in addition, don't depend on knowing an accurate source flux or a noise diode calibration standard 
(cal) temperature. 

The robust parameters specify only the beam shape using ratios of the first sidelobe and 
main beam. The peak power in the first sidelobe, as a fraction of the peak of the main beam, 
is Pn,FS,u,r-, the corresponding quantity for the main beam, Pn,MB,u,i is unity by definition. We 
obtain 



Pn,FS,u,i = 0.0175 (15a) 

is 



and the ratio of beam efficiencies as defined by equation 13 

VFS \ 



VMB 



)u. 



0.0861 , (15b) 



Obviously, these two parameters are not independent because, to some extent, they measure the 
same thing. 

In the same way we can compute three non-robust parameters, 

Half-power beam width: 

{HPBW)ui = 0.597 ^^ arcmin. (16a) 

KperJy^i^ 

Main beam efficiency: 

{mB)u,i = 0.840 , (16b) 



And, combining equations 16b| and 15b| to get the total efficiency in the main beam and first 
sidelobe: 

{ms + VFs)ui = 0.840(1 + 0.0861) = 0.912 . (16c) 



To calculate HPBW in equation 16a we use diameter d^s = y ^"^ and the HPBW for a 
uniformly-illuminated aperture, HPBW = l.OS-j — . For the data, the observed HPBW is derived 
from the Gaussian fit, and the value expected from the model is derived from the effective area 
AeEi which is in turn derived directly from KperJy, namely from observations of calibration sources 
with known flux using a standard noise diode. This dependence on KperJy is what makes these 
parameters non-robust, because KperJy is affected by telescope imperfections (surface roughness 
and subreflector misalignment), loss in waveguide feeds, errors in source calibrator flux, and errors 
in noise diode calibration temperature. 

Our beam parameterization is done through Gaussian fits to the data, and thus we should 
compare those fits with Gaussian fits to the ideal uniform illumination pattern. Gaussian fits don't 
represent the theoretical model perfectly, but we can compute correction factors Pps^ Eps^ and 
Emb, and H for the Gaussian fits to the ideal uniform circular aperture. As we described above, 
the Gaussian-fit derived parameters are subscripted only with "[/" . These correction factors are 



'^™"'" (17a) 



H = 


HPBWu,i 




'nMB,U 




rjMB^U,! 


Eps — 


VFS,U 
VFS,U,I 


Pt:^^ — 


( Pn,FS,U 



p. 



n,FS,U,I 



(17b) 

(17c) 

(17d) 



For example, multiplication of the ideal HPBWjj,! by H gives HPBWu, the value that is derived 



from a Gaussian fit to the ideal beam shape; the factor 0.597 in equation 16a is multiplied by 
0.961, becoming 0.574. The ideal beam is wider at the top and falls more rapidly than the 
Gaussian; this makes H < 1 for a uniformly-illuminated aperture with no blocking. 



Below, in |4.3| , we will consider the effects of blockage. There are corresponding correction 
factors for blockage, which we quote here for convenience. Table ^ gives values of the {H, P, E) 
factors in equation ^ for the unblocked case (0%) and for two different blocked cases with 10% 
and 20% blockage. Empirically, the factors change slowly with blockage, so that these few tabular 
entries are sufficient. 



4.2. Comparison of Observations with the Unblocked Aperture 

Table |3| compares some of the observed parameters with those predicted for an unblocked 
aperture (equations 15 and |l^). Arecibo's beam doesn't agree at all with these values. In 



particular, consider the two robust parameters Pn,FS,u and ( ^^^ j . The observed values are 



higher than the unblocked values by factors ~ 2.2 and 4.6, respectively. Conventional feeds, 
including those used here, produce a tapered illumination, which decreases the sidelobe level 
relative to the uniformly illuminated case. It is impossible to increase the sidelobe level without 
introducing a relative decrease in the center of the illumination pattern — a donut-like pattern. 
For Arecibo's Gregorian, this occurs naturally because of blockage of the primary reflector by the 
central support structure. We can use Table ^ to estimate the fraction of blocked area for LBW. 



4.3. The Effects of Blockage 

We calculate blockage by subtracting the voltage pattern of a circular blocking aperture from 
that of the illuminated aperture, taking care to keep the power ratio equal to the area ratio. We 
keep the effective area constant, because A^s is what our observations provide. We normalize the 
voltage at the center of the effective area's pattern to unity and plot the above quantities versus 
effective area ™ Figure |2|. For comparing with observations we assume that the blocked radiation is 
scattered into the sky, not absorbed, so that equation ^ holds; this reduces tjmb by ni^^^^ted'^^^en 
relative to the assumption that the blocked radiation is absorbed. Recall that, according to 
Babinet's principle, the effect of geometrical blockage is doubled; below we will derive blockage 
^ 10%, which corresponds to physical blockage ~ 5%. 

Figure ^ exhibits plots of the parameters versus the fractional blocked area. Xs on these 
plots show the measured quantities; the points are placed vertically using the observed value 
and horizontally using the curves. Thus, the "derived" blockages effective area ^^^* correspond to 
the model can be read off of the horizontal axes. For the robust parameters, the blockages lie 
somewhat below and above 0.10 in the top two frames, respectively. It seems reasonable to adopt 
a blockage of 0.10 for discussion purposes. 

For this particular blockage. Figure y exhibits the beam shape and, also, the Gaussian fit. The 
Gaussian is a reasonably accurate rendition of the theoretical model, but has its faults: it doesn't 
fall to zero between the main beam and the sidelobe and it doesn't reproduce the details of the 
shape very well. Then again, neither the Gaussian nor the model fits the data particularly well. 
In fact, they cannot possibly do so because the sidelobe structure changes with (p, while the model 
is circularly symmetric. An additional limitation of both our parameterization and our model is 
the neglect of additional sidelobes, which contribute non-trivially to the total beam efficiency. 

This model of a uniformly illuminated blocked aperture is useful as an indication of expected 
properties of the main beam and sidelobes. However, it is limited and inaccurate for the following 
reasons: 

(1) Most feed systems don't provide uniform illumination; rather, the illumination is tapered. 
Tapering increases HPBW and r]MB and decreases sidelobe levels, which is opposite to blockage. 
Thus, the actual geometrical blockage must be somewhat larger than our adopted 10%. Because 
these two effects work in opposite directions, we surmise that the relationships among our derived 
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parameters (such as HPBW and rj) remain approximately correct. 

(2) As a telescope is pushed to higher frequencies, imperfections in the telescope surface 
become increasingly important. If these imperfections are random, then they scatter radiation 
away from the main beam according to the Ruze (1966) formulation and will decrease KperJy and 



all beam efficiencies rj. Using this reduced KperJy may make equation 16a provide a HPBW that 
is too large. The ratio of FS to MB should not be affected much unless the correlation lengths of 
the surface imperfections are large — but, of course, for most telescopes, there do exist some large 
scale correlations. 

(3) At higher frequencies, large-scale imperfections in the telescope surface and, also, in 
the mirror alignment become important for most telescopes. These produce significant changes 
in near-in sidelobes. These effects decrease KperJy, increase HPBW, tend to increase the first 
sidelobe power, and tend to transfer some of the main beam power to near-in sidelobes. We expect 
Vmb + Vfs to decrease, but their power to go primarily into near-in sidelobes so that the antenna 
temperature for an extended source will not decrease very much. Of course, the detailed behavior 
is complicated, and in addition small-scale random surface irregularities also become important. 



4.4. A More Accurate Determination of Blockage 



Equation 16a describes a non-robust relationship between two robust parameters, HPBW and 
A; these parameters are robust for frequencies below those where large-scale surface irregularities 
affect the main beamwidth. This equation is non-robust because it depends on KperJy. This 
particular value of KperJy is that for a telescope with perfectly smooth reflecting surfaces and 
perfect mirror alignment. This value can be calculated by the telescope and feed engineers; we 
denote it by KperJy i^eai- It can also be measured at frequencies low enough for the telescope to be 
regarded as perfect; however, such measurements rely on knowing accurate source fluxes and cal 
temperatures. 



We can observationally determine the blockage by first rewriting equation ICa: 



(0.597/?) 



HPBWa 



(18) 



Here /9 is a parameter that accounts for the blockage and is given hy (3 = H (upg^) , where H 
comes from Table |2| and the HPBW ratio from Figure pi. Next we measure HPBW; knowing 
KperJy ideal, we obtain the blockage factor (3. 

In practice, one measures HPBW over a range of wavelengths, using feeds with identical 
illumination patterns, least-squares fit for — '■ — j4-, and calculates /3. We obtained results at 

Arecibo for A = 6 ^ 70 cm; the fit is very good and, assuming KperJyi^eai = 10.5, it yields 
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P = 0.88. From Figure |2|, this corresponds to about 13% blockage. We believe that the blockage 
so derived should be accurate, but it does depend on knowing a reliable value for Kp&-Jy ideal- 



5. STRUCTURE IN THE FIRST SIDELOBE 



We introduced the Fourier decomposition of the first sidelobe in section p.3| ; we obtained a 
series with 8 points and we use iV = 6 to reconstruct the sidelobe. This allows us to reconstruct 
four Fourier components, numbered zero to three. Figure Q exhibits the az dependence of these 
components, and we briefly discuss these data to illustrate the power of this representation, 
because the {az, za) dependences of these Fourier coefficients reveal their probable production 
mechanisms. 

The average amplitude (the zeroth Fourier component, top panel) seems to increase slightly 
near transit, which probably reflects the slightly increased blockage that occurs there. 

The angles of maximum response (j)max of the first and third Fourier components exhibit a 
striking dependence with az: the points follow the slope of the dotted line, meaning that these 
components are fixed with respect to the ground. The third Fourier component would respond 
to blockage by the triangle structure, which is a large, major structural entity that supports the 
feed structure; its az dependence implies that it is doing just that. The interpretation for the first 
Fourier component is less obvious. 

The position angles of the second Fourier component are independent of az, implying that 
they are fixed with respect to the azimuth arm. This is consistent with the blockage of the azimuth 
arm, which is a linear structure roughly centered on the illumination pattern and should produce 
a second Fourier component. 

The first and third Fourier components are fixed with respect to the ground, in contrast to 
the behavior of the coma lobe, which is fixed with respect to the azimuth arm. We conclude that 
the sidelobe structure and coma lobe are, for the most part, unrelated. We reaffirm this conclusion 



in §3.3, where we discuss the sidelobe polarization. 



6. THE POLARIZED STOKES PARAMETERS 

One derives the beam structure of the polarized Stokes parameters by observing an 
unpolarized source, or by subtracting off the polarization of a polarized source to simulate an 
unpolarized source. An unpolarized source contains power in all polarizations, so beam structure 
in the polarized Stokes parameters reflects enhanced gain of one polarization over its orthogonal 
counterpart. Thus, in this section we are discussing the response of the telescope in Stokes Q, U , 
and V to an unpolarized source. Note that this response, unlike the Stokes / response, can be 
positive or negative. We assume that the on-axis polarization has been corrected as described in 
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the companion paper Heiles et al (2001), so that (Q, U,V) = at beam center. 



6.1. Squint and Squash: Parameters for the Polarized Main Beam 

For each polarized Stokes parameter such as Q, we express the response of the main beam 
with the first two terms of a two-dimensional Taylor expansion. The first derivative term, 
conventionally called "beam squint" , is odd and corresponds to displacement of the beam center 
in the two orthogonal polarizations. Beam squint is characterized by an angular magnitude and a 
direction and has a two-lobed structure on the sky, with one positive and one negative lobe. The 
second term, which we call "beam squash" , is even and corresponds to a difference in beam width 
between the two polarizations. Beam squash is characterized by a magnitude and an orientation 
and has a four-lobed structure, with lobes on opposite side of beam center having identical signs. 
We assume that the coma lobe is unpolarized. 

To be more precise, consider Stokes Q as an example, which is polarization X minus Y. We 
define the normalized power pattern for a polarization X as ^ ^ .q| Then the power received from 
the source in each polarization is 



Px = Sx exp 



Py = Sy exp 



-1-2^1^ "coma Oo 



e+^y 



u 2 J y'^ "coma 00 J 



(e-^^' 






(19a) 



(19b) 



where Sx and Sy are the source fluxes in the two polarizations. Recall that Q, 6q, and Qq all 
depend on (p. The unsubscripted P is the sum (Stokes /), while the difference gives the measured 
Stokes Q: 



Pc 



QPn + \ 



dPn 

den 



dP 



(20) 



The quantities 56q{(J)) and 6QQ{(f>) are not themselves the beam squint and squash. Rather, 
the squint and squash are represented by (magnitude, position angle), {0 squint i4> squint) and 
{0 squash A squash), and are given by 



"squint COS\ 



^squint) 



(21a) 



Recall that above we defined the unsubscripted G (which includes Gmax) to apply to 
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SQq = Qsquash COS 2(0 - (psquash) ■ (21b) 



The quantity 59q varies as cos (j) because it is like a pointing offset and has a direction. In contrast, 
the quantity 6Qq varies as cos 2(p: in beam squash, the beamwidth varies as cos 20 so it only has 
an orientation. 



6.2. Results for Beam Squint and Squash 



We least-squares fit equation 20 to the data, regarding P„ as known. This is acceptable, 
even though Pn is itself a result of a least-squares fit, because its data are Stokes / with high 
S/N. We do this in all three polarized Stokes parameters. Figures |5| and |^ show representative 
plots for Stokes {Q,U). The squint and squash amplitudes are in the few arcsecond range and 
there is significant variation, but no easily interpretable physical dependence, on (za, az). On the 
other hand, Stokes V (Figure j^) is more smoothly behaved: the squint is somewhat smaller and 
maintains a constant 4>, while the squash is very small. 

For the linearly polarized Stokes parameters {Q,U), we expect beam squash but no squint 
(Tinbergen 1996). The prediction regarding squints is not satisfied at all. Instead, beam squints for 
Stokes {Q, U) are comparable to that for V . The behaviors of (psquint for {Q, U) are uninterpretable 
in simple terms and provide no hint about the cause of the squint. The beam squash amplitude 
variations are also uninterpretable. The angles (jisquash should be fixed with respect to the azimuth 
arm; very roughly speaking, they are. However, they should differ by 45° between Q and U; they 
don't. 

For Stokes V, one expects beam squint but no beam squash (see discussion by Troland 
and Heiles 1982). This prediction is reasonably well satisfied. The squint amplitude is roughly 
constant and its angle (j)squint is constant, fixed with respect to the azimuth arm as it should be. 
The squash is small, far smaller than for Stokes (Q, [/), although it is not zero. 

The reasons for the severe departure of Stokes (Q, U) from theory are completely unclear to 
us. The absence of a clearly defined dependence with on az suggests a combination of effects that 
are fixed with respect to the ground and the az arm; the obvious candidates are aperture blockage 
and secondary/tertiary mirror misalignment. 



6.3. Fourier Parameters for the Polarized Sidelobes 

For each polarized Stokes parameter such as Q, we follow §^ and express the response of the 



sidelobe in terms of Fourier components. Figures |g, y, and 11 display these coefficients for the 
three polarized Stokes parameters. The zeroth Fourier coefficient describes the net polarization 
averaged over the sidelobe ring and is shown in the top panels. The first sidelobe is strongly 
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linearly polarized. 

In Stokes Q and U, the first Fourier component dominates and sometimes produces fractional 
polarizations — — j-^ — ^ 0.5. Their angles of maximun response (f)max have a complicated 
behavior: their slopes with respect to az indicate that they are approximately fixed with respect 
to the ground, but (jjmax suffers a 180° jump near az = (—20°, 20°) for {Q,U) respectively. The 
dependences for the second and third Fourier components are even less clear. 

Because of the dominance of the first Fourier component, the sidelobe linear polarization 
has a strong component that is fixed with respect to the ground (apart from the 180° jumps), in 
contrast to the behavior of the coma lobe, which is fixed with respect to the azimuth arm. Thus 
we reaffirm the conclusion of §^ the sidelobe structure and coma lobe are, for the most part, 
unrelated. 

The angles shown in Figures |8|, P, and 11 are not the position angles of linear polarization. 
Rather, they are the angles (j), defined in Figure |^, at which the particular Fourier coefficient 



is maximum. In contrast. Figure 10 exhibits the position angle of linear polarization 
P^iinpoi = 0.5tan^^^ for the four Fourier coefficients. As expected, PAn^poi for the zeroth 
component follows the azimuth arm in its rotation with respect to the sky. For the other 
components the behavior is more complicated, which occurs because 4>max changes with azimuth: 
in other words, the (j) of maximum response moves around the sidelobe ring differently in Q and U 
as the azimuth changes, which leads to rotation of PA^npoi with respect to the azimuth arm as it 
rotates when tracking a source. 

In Stokes V, the first Fourier component dominates and produces fractional circular 
polarization -j ~ 0.1. For the first and third Fourier components the angles of maximun response 
(pmax are nearly constant, being fixed with respect to the ground. The dependence for the second 
Fourier components is unclear, perhaps because it is so weak. 



6.4. The Effects of Beam Squint and Squash on Extended Sources 

Suppose one is observing a large-scale feature where the brightness temperature Tg varies 
with position. One can express this variation by a two-dimensional Taylor expansion. The first 
derivative is described by a single term with a magnitude and direction; the second derivative 

^2^7^ /7'^T^ rf^T' 

is described by three terms, ^^J! , -j^, dazdza ' ^ach of which is an ordinary second derivative 
with a direction. In this section and its figures, Tb refers to Stokes / so is t^vice its 
conventionally used value. 

Beam squint, by its nature, responds to the first derivative and only slightly to the second; 
beam squash responds primarily to the second. To illustrate the effects of polarized beam 
structure, we present plots of the response of the polarized beam to these derivatives. In these 
plots we adopt derivatives of 1 K arcmin"^ and 1 K arcmin"^ to facilitate scaling to other 
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measured derivatives, not because these values are necessarily realistic; also, in these figures the 
brightness temperatures refer to Stokes /, not 2, so they are twice the conventionally used values. 
The response comes from three separable causes: main beam squint, main beam squash, and the 
first sidelobe polarization structure. We show these responses individually and, also, we show the 
total response. 



Figures 12 and O show Stokes Q and U . In these two figures, the position angles PA are 
defined on the sky with the conventional astronomical definition. For the first-derivative case, the 
dominant contribution comes from the first sidelobe, which contributes far more than the main 
beam squint. The amplitudes are roughly constant with az. The angles of maximum response 
are, very roughly, fixed on the sky (they change less rapidly than az), which is consistent with the 
first sidelobe's effects being caused by blockage. For the second derivative case, the first sidelobe 
contribution dominates over main beam squash, but not by much. 

Figure ^ shows the Stokes V result. In this figure, PA is the angle of maximum response 
defined relative to the azimuth arm. The total amplitudes for both derivatives are smaller by 
a factor ~ 10 than for Stokes {Q,U). For the first derivative, the PA is independent of az, 
as expected from the squint result. The sidelobe is again the dominant contribution, but not 
overwhelmingly so: in fact, its amplitude is about twice that of the main beam squint but the 
direction is opposite. This has the curious effect of making the amplitude of the total beam 
response about the same with and without the sidelobe — but reversing the direction of gradient for 
maximum response! For the second derivative, the first sidelobe and main beam squash contribute 
comparably and the PA behavior is correspondingly complicated. 

We emphasize that these particular results can produce serious systematic errors when 
measuring the polarization of extended emission. The results represented in Figures 12, ^, ^ 



are contributions in the polarized Stokes parameters {Q,U,V) that arise from spatial gradients 
in the total intensity Stokes parameter /. The fractional polarization of extended emission tends 
to be small, so spatial gradients in / can be very serious. Correcting for them at Arecibo is a 
complicated business because of the PA variation with {az,za). It is also an uncertain business, 
especially for {Q,U) and somewhat less so for V, because these variations are unpredictable and 
must be determined empirically. Presumably, corrections at conventional telescopes are more 
straightforward. 

In that spirit, these plots are more useful as guides of what can be reliably measured at 
Arecibo than they are for correction of the measurements. Consider, for example, measuring 
Zeeman splitting of the 21-cm line in emission, which involves measuring Stokes V. Figure 14 shows 
that a brightness temperature gradient -^ = 1 K arcmin^^ produces an artificial temperature 
Vart ~ 0.04 K response in Stokes V; let this ratio be denoted by /C = ^y Lg ~ ^-^^ arcmin. Then 
it's easy to show that, for a Lande g factor 2.8 Hz /iG~^, the corresponding artifical magnetic field 
is Bart = 1690/C^ /xG, where -^ is the velocity gradient in km s~^ arcmin"^ and /C is in arcmin. 
For a representative 3I = 1 km s~^ deg~^ we get S^rt ~ 1-1 /^G. Heiles (1996) presents statistics 
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of the velocity gradient measured with a 36 arcniin HPBW; gradients of 1 km s~^ deg~^ are not 
uncommon. Gradients might be larger when measured with smaller HPBW. Typical values of B 
are in the fiG range, so this effect can be — but is perhaps not always — serious! 



7. IMAGES OF THE BEAM IN ALL STOKES PARAMETERS 

Expressing the beam properties in terms of parameters is especially appropriate for 
quantitative work. However, for getting an intuitive feel for the situation there's nothing like an 
image. So here we provide images for one of our several hundred beam maps. As with our previous 
detailed figures, these images are all for LEW (1175 MHz) observing B1749+096; the source was 
setting near Arecibo's za limit of 20°, with (az, za) ~ (64°, 18.4°). 

The main beam parameters are as follows: 

For Stokes /, (60,61) = (4.00,0.36) arcmin (both are HPBW); 0beam = 91.1°; acoma = 0.048; 
0coma = 41.4°; mean height of first sidelobe= 0.029. 

For Stokes Q, {0 squint, 4) squint) = (0.019', -173°) and {Q squash A squash) = (O.ll', 19°). As in 

Figure |5|, Q is oriented 45° to the azimuth arm. 

For Stokes U, {9 squint, (k squint) = (0.056', -42°) and {Q squash, (t> squash) = (0.061', -47°). As in 

Figure |6|, U is aligned with the azimuth arm. 

For Stokes V, {0 squint, (k squint) = (0.045', -10°) and {Q squash, (t> squash) = (0.007', 37°). 

The beam is considerably elliptical, which is a result of the high zenith angle. The Stokes 
Q beam shows the classic four-lobed pattern, beam squash, for linear polarization; in contrast, 
the U beam is polluted with a large unpredicted squint and a huge linearly polarized sidelobe. 
Theoretically, the beam squashes in (Q, U) must be aligned with their orientations, and in 
particular the two four-lobed patterns should be rotated 45° with respect to each other; this is not 
the case! The V beam is reasonably well-behaved. We commented on these aspects in more detail 
above. 



It is a pleasure to acknowledge insightful and helpful comments by the referee, which 
substantially improved several aspects of this paper. This work was supported in part by NSF 
grant 95-30590 to CH. 
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8. FIGURE CAPTIONS 

Fig. 1. — Grey-scale/contour image of a main beam having Ocoma = 0.2, i;^coma = 22.5°, Qq = 3.4 
arcmin, 0i = 1 arcmin (both 0's are HPBW in this figure), and i;^beam = —67.5°. The dashed 
lines intersect the true beam center; the dashed-dot hnes intersect the assumed (from the telescope 
pointing) beam center, and are the paths of the telescope in the observing pattern. Solid contours 
are (0.1, 0.2, . . .) of the peak; dashed contours are (0.01,0.02, . . .). 

Fig. 2. — Blockage-sensitive, measurable parameters versus the ratio of blocked area to effective 
area. Solid lines show the parameters as derived by fitting the main beam and first sidelobe to 
Gaussians; dotted lines show the actual parameters for the theoretical model. The top two frames 



show the robust parameters of equation 15 and the third and fourth show the non-robust parameters 



of equation |16 



Fig. 3. — Normalized power pattern P„, and also the integrand OP^, for (1) the standard model of a 
uniformly illuminated blocked aperture (solid, dash-dot lines), and (2) its Gaussian-fit counterpart 
(dash, dash-dot-dot-dot lines). The squares and diamonds are representative data points for P„ 
from the LBW feed at 1415 MHz, obtained by averaging different cuts in one single observing 
pattern. 

Fig. 4. — az dependence of the first four Fourier components of the first sidelobe in Stokes / for 
LBW at 1175 MHz. The dotted lines correspond to a component being fixed with respect to the 
ground. Solid lines are amplitude and the points are position angle 4>max-, which is the angle (/> at 
which the component peaks {4> is defined in Figure |l|). For the nth Fourier component (n = ^ 3), 
each point occurs n times over each 360° interval in 0, and all of those possibilities are plotted to 
better clarify the az dependence. 

Fig. 5. — {za, az) dependence of amplitudes (in arcmin) and {(j) squint -.(k squash) of beam squint and 
squash for Stokes Q. (Here (j) is labeled "PA"). The dotted line corresponds to being fixed with 
respect to the sky or ground. In this figure, Q is oriented 45° to the azimuth arm and 
fixed w^ith respect to it. 

Fig. 6. — (za, az) dependence of amplitudes (in arcmin) and [(j) squint -.(^ squash) of beam squint and 
squash for Stokes U . (Here (j) is labeled "PA"). The dotted line corresponds to being fixed with 
respect to the sky or ground. In this figure, U is aligned with the azimuth arm and fixed 
w^ith respect to it. 

Fig. 7. — {za, az) dependence of amplitudes (in arcmin) and {(f) squint ■,<i> squash) of beam squint and 
squash for Stokes V. (Here (j) is labeled "PA"). The dotted line corresponds to being fixed with 
respect to the sky or ground. 

Fig. 8. — az dependence of the first four Fourier components of the first sidelobe in Stokes Q for 
LBW at 1175 MHz. The dotted lines correspond to a component being fixed with respect to the 
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ground. Solid lines are amplitude and the points are position angle (pmax, which is the angle (j) at 
which the component peaks ((/> is defined in Figure |l|). For the nth Fourier component (n = ^ 3), 
each point occurs n times over each 360° interval in (j), and all of those possibilities are plotted to 
better clarify the az dependence. 

Fig. 9. — az dependence of the first four Fourier components of the first sidelobe in Stokes U for 
LBW at 1175 MHz. The dotted lines correspond to a component being fixed with respect to the 
ground. Solid lines are amplitude and the points are position angle (pmax, which is the angle (j) at 
which the component peaks ((/> is defined in Figure |l|). For the nth Fourier component (n = ^ 3), 
each point occurs n times over each 360° interval in (j), and all of those possibilities are plotted to 
better clarify the az dependence. 

Fig. 10. — az dependence of the position angle of linear polarization PAnnpoi for the first four 
Fourier components of the first sidelobe for LBW at 1175 MHz {PA is the conventional astronomical 
definition). The dotted lines correspond approximately to PAnnpoi being fixed with respect to the 
azimuth arm. 

Fig. 11. — az dependence of the first four Fourier components of the first sidelobe in Stokes V for 
LBW at 1175 MHz. The dotted lines correspond to a component being fixed with respect to the 
ground. Solid lines are amplitude and the points are position angle (pmax, which is the angle cp at 
which the component peaks {(p is defined in Figure |l|). For the nth Fourier component (n = ^ 3), 
each point occurs n times over each 360° interval in (j), and all of those possibilities are plotted to 
better clarify the az dependence. 

Fig. 12. — az dependence of the telescope's response to the first spatial derivative in brightness 
temperature (top two panels) and to the second derivative (bottom two panels), all for Stokes Q. 
PA is direction of maximum response. In this figure, Q and PA are defined on the sky with 
the conventional astronomical definition. Tg refers to Stokes I so is twice is conventionally 
used value. 

Fig. 13. — az dependence of the telescope's response to the first spatial derivative in brightness 
temperature (top two panels) and to the second derivative (bottom two panels), all for Stokes U. 
PA is direction of maximum response. In this figure, U and PA are defined on the sky with 
the conventional astronomical definition. Tg refers to Stokes I so is twice is conventionally 
used value. 

Fig. 14. — az dependence of the telescope's response to the first spatial derivative in brightness 
temperature (top two panels) and to the second derivative (bottom two panels), all for Stokes V. 
PA is the (j) of maximum response, defined relative to the azimuth arm as in Figure |l|. Tb 
refers to Stokes / so is twice is conventionally used value. 

Fig. 15. — Grey-scale/contour of the Stokes I main beam and first sidelobe reconstructed 
from the derived parameters. Solid contours are (0.1,0.2,...) of the peak; dashed contours are 
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(0.01,0.02,...). 

Fig. 16. — Grey-scale/contour of the Stokes Q main beam and first sidelobe reconstructed from 
the derived parameters. Black contours are for areas with negative Q with the grey scale tending 
towards white; white contours are positive Q with the greyscale tending towards black. Contours 
are in percent of Stokes / at beam center and spaced by 0.4%; the 0% contour is omitted. As in 
Figure |5|, Q is oriented 45° to the azimuth arm. 

Fig. 17. — Grey-scale/contour of the Stokes U main beam and first sidelobe reconstructed from 
the derived parameters. Black contours are for areas with negative U with the grey scale tending 
towards white; white contours are positive U with the greyscale tending towards black. Contours 
are in percent of Stokes / at beam center and spaced by 0.4%; the 0% contour is omitted. As in 
Figure |6|, U is aligned with the azimuth arm. 

Fig. 18. — Grey-scale/contour of the Stokes V main beam and first sidelobe reconstructed from 
the derived parameters. Black contours are for areas with negative V with the grey scale tending 
towards white; white contours are positive V with the greyscale tending towards black. Contours 
are in percent of Stokes / at beam center and spaced by 0.2%; the 0% contour is omitted. 
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Table 1. REPRESENTATIVE OBSERVED PARAMETER VALUES 



RCVR FREQ 



SOURCE 



P^,i 



vfs 

VMB 



KperJy 



HPBW 



rjMB 



rjMB + Vfs 



430G 


430 


B1634+269 


10 


.039 


0.33 


10.3 


10.9 


0.66 


0.88 


LBW 


1175 


B1634+269 


10 


.043 


0.32 


8.7 


4.0 


0.56 


0.75 


LBW 


1415 


B1634+269 


10 


.046 


0.32 


7.5 


3.4 


0.50 


0.67 


LBW 


1666 


B1634+269 


10 


.042 


0.31 


6.9 


2.9 


0.48 


0.62 



Note. 



za is in degrees and HPBW is in arcmin. 
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Table 2. FACTORS IN EQUATIONS |T^ VERSUS 
BLOCKAGE 



FACTOR 



0% 



10% 



20% 



H 


0.961 


0.955 


0.952 


Pfs 


1.038 


1.038 


1.033 


Emb 


1.060 


1.067 


1.070 


Efs 


0.865 


0.925 


0.947 
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Table 3. COMPARISON OF OBSERVED AND UNBLOCKED MODEL PARAMETERS 



RCVR 


FREQ 


Pfs 
(Pfs)u 


(ifs/imb) 


dcB 


HPBW 


nuB 


VMB+VFS 


(vfs/vmb)u 


(HPBW)u 


iVMB)u 


(VMB+VFS)U 


430G 


430 


2.2 


4.7 


190 


0.87 


0.74 


0.92 


LBW 


1175 


2.4 


4.6 


175 


0.80 


0.63 


0.79 


LBW 


1415 


2.5 


4.6 


162 


0.75 


0.56 


0.70 


LBW 


1666 


2.3 


4.4 


156 


0.71 


0.54 


0.65 



Note. — See « for parameter definitions, dcff is in meters. 
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SIDELOBE HEIGHT INCREASES WITH BLOCKAGE 
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10% BLOCKED APERTURE, ITS GAUSSIAN EIT, AND SOME DATA 

SOLID, THEORETICAL MODEL BEAM RESPONSE 
DASHED, GAUSSIAN EIT BEAM RESPONSE 
DASH-DOT, THEORETICAL MODEL BEAM INTECRAND 
DASH-DOT-DOT-DOT, GAUSSIAN EIT BEAM INTEGRAND 
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bw 1175MHz BRDO B1749+096 1 2-SEP-2000 STOKES Q 
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